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Giesbertz, Gritsenko and Baerends (GGB) have stated 
that the occupation numbers nk in time-dependent den- 
sity matrix functional theory are time independent in 
the "adiabatic" approximation (AA) for any ground- 
state functional [l[ (see also Q). It is important to 
know whether this statement is true as it has implica- 
tions for the design of functionals capable of generating 
time-dependent occupation numbers. Here we show that 
the argument given by GGB to support this statement 
is incorrect. The statement, however, is true; it follows 
uite generally from the stationarity of the ground state 



The equation of motion for the one-body reduced den- 
sity matrix 7 implies idnk /dt = W], k — Wkk [2| , where 



Wki = y^^w kqrs T 
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with T srq i = {^IcjcJcrCsl^} and 



Wkqr 



dxdx'cj>l (x)<f> q ( x ')——fa(x')fa (x) . (2) 



In the AA, the memory-dependent functional T([7]; t) on 
the right-hand side of Eq. (JT]) is approximated by the 
ground-state functional To [7] evaluated for 7(i). GGB 
argue that the invariance of the ground-state interaction 
energy functional Wo — Wo [7] with respect to the change 
fa, —> e lak <j)k in the phases of the natural orbitals implies 
dnk/dt — 0. Therefore, they claim to prove the implica- 
tion dWo[y]/dak = =$> dnk/dt = 0. The crux of their 
argument is the statement 
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where Wo.kk are defined in the same way as the Wkk but 
with ground-state quantities. To establish Eq. ©, GGB 
use the identity 
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da k J H* k {x) 
and the statement 
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quoted from Ref. 
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where it was derived from 
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Here W [y] = |niin {€?} £ fe(?rs ^kqrs^srqk^p) and 
parametrize a constrained search over A^-representable 
T srq k that contract to 7. It was argued 0] that Eq. (O 
follows from Eq. © because the first term vanishes due 
to the variational nature of the constrained search. But 
there are two flaws with this argument: (i) Eq. ^ it- 
self is manifestly incorrect because Wo[-f] has no {£ p } 
dependence after the constrained search has been per- 
formed: the operations ^2 p and min{^} do not com- 
mute; (ii) The variational character of W)([0i];£ p ) = 
\ Hkqrs Wkqrs? srqk{£, P ) at the minimizing £ P for fixed 
{fa,rii} does not imply that the gradient with respect 
to £ p is zero, because Wo([</>i]; £ p ) is only stationary with 
respect to the subspace of degrees of freedom that 
are orthogonal to the fa degrees of freedom as the latter 
are constrained. Ultimately, the argument is incorrect be- 
cause it does not account for the fa dependence of T srq k . 

For the specific case of approximate Wo[ni, fa, </>*] that 
contain only Wkqkq and Wkqqk Coulomb integrals and in 
which the T srq k are functions of Hi, it might seem that 
Eq. ([5]) can be verified by an explicit calculation of the 
functional derivative. However, such a calculation is not 
valid because the variation (f)* k — > <j)* k + 8cf)* k holding fixed 
all other fa* and all cf>i corresponds to a non-Hermitian 
7 + Sj. Hence, such a variation goes outside the physical 
domain of Wo [n*, fa, <j>*\. The functional derivative of 
Wo [7] with respect to an orbital should be understood as 
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dxdx'fa [x] — -j—p^-rfa {x') , 
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where now Wo [7] = 5/ dx\dx 2 ^^ Tq{[^]\x\X2,x 1 X2) 
and To{['y\,XiX2,x l 1 x l 2) is the ground-state two-body re- 
duced density matrix functional jB[. Clearly, Eqs. ((H) and 
([7]) cannot justify Eq. (J3)) because the right-hand side of 
Eq. (j3|) depends on degrees of freedom oiTo{xiX2,x' 1 x' 2 ) 
that are integrated out in the definition of Wo [7]. This 
information cannot be recovered by taking the derivative 
with respect to oik- 
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